clear all
clc;
[x,y] = runge(@dyfun, [0 8.032], [9.279 7.2 0 7.2 0 0.27 11.076], 0.001);
figure(1)
plot(y(5,:),y(6,:));
title('弹道曲线')
xlabel('X/m');
ylabel('Y/m');
grid on
figure(2)
plot(x,y(1,:))
title('速度曲线')
xlabel('t/s');
ylabel('V/(m/s)');
grid on
figure(3)
plot(x,y(2,:))
title('弹道倾角曲线')
xlabel('t/s');
ylabel('bt/°');
grid on
figure(4)
plot(x,y(4,:)-y(2,:))
title('攻角曲线')
xlabel('t/s');
ylabel('alpha/°');
grid on
figure(5)
plot(x,y(6,:))
title('时间高度曲线')
xlabel('t/s');
ylabel('Y/m');
grid on